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Abstract 



> 

(yQ ' Recent angle resolved photoemission (ARPES) results for the insulating 

o ' 

cuprate S^CuC^Ck) have provided the first experimental data which can be 

O " 
in 



directly compared to the (theoretically) well-studied problem of a single hole 



propagating in an antiferromagnet. The ARPES results reported a small 

a" 

bandwidth, providing evidence for the existence of strong correlations in the 



cuprates. However, in the same experiment some discrepancies with the famil- 
iar 2D t — J model were also observed. Here we discuss a comparison between 



^ ■ the ARPES results and the quasiparticle dispersion of both (i) the t — t' — J 



Hamiltonian and (ii) the three-band Hubbard model in the strong-coupling 
limit. Both model Hamiltonians show that the experimentally observed one- 
hole band structure can be approximately reproduced using reasonable values 
for t', or the direct oxygen hopping amplitude t pp . 

To appear in Phys. Rev. B, April 1st, as a Rapid Communication. 



1 



Typeset using REVTeX 



The dispersion of the hole quasiparticles in the normal state of the two-dimensional 
(2D) Cu02 planes is of crucial importance in the search for a microscopic theory of the high 
temperature superconductors. To be specific, comparison between experimentally observed 
and theoretically predicted quasiparticle dispersions allows for the scrutiny of Hamiltonians 
proposed to describe the carriers in these planes. Recently, angle resolved photoemission 
(ARPES) data have shown that for Bi2212 and YBCO, both near optimal doping, the 
quasiparticle bandwidth is only ~ 0.3eV. In addition, the dispersion contains an interesting 
flat region around X = (n, 0) and Y = (0, 7r) (in the notation of the 2D square lattice). The 
apparent universality of these results among different hole-doped compounds implies that 
their origin must lie in the CuC>2 planes. The experimentally observed small bandwidth sug- 
gests that strong correlations are crucial for a proper description of carriers in the cuprates. 
Further, calculations based upon the t — J model || have shown that the flat regions near 
X and Y are naturally obtained for holes moving in an anti ferromagnetic (AFM) back- 
ground, the latter being a direct consequence of the effect of strong correlations at half 
filling. Then, for doping levels such that the AFM correlation length, £af, is robust, the 
dispersion obtained from the 2D t — J model should approximate the observed quasiparticle 
dispersion. 

In a recent letter, Wells et al. |§ reported ARPES measurements on S^CuC^C^. This 
material is an insulating layered copper oxide that is difficult to dope. These results are 
extremely important because this experiment precisely corresponds to the well-studied theo- 
retical problem of a single hole moving in a quantum antiferromagnet, and thus is a problem 
that serves as a test case for any proposed microscopic Hamiltonian, as we now discuss. 
From Ref. || the main message for theorists can be summarized in two items: (i) on the one 
hand, the ARPES spectra shows that the 2D t — J model accurately describes the results 
along the diagonal (0, 0) to (tt, it) in momentum space. The state with the highest energy 
in the valence band is obtained at M = (7r/2,7r/2) as expected from many calculations [|J. 
Further, the bandwidth of the quasiparticle is very close to the bandwidth obtained for one 
hole in the 2D t — J model using various techniques [||||. This result is of considerable 
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importance for theories of high-T c as it supports the use of one-band models to describe 
the electronic properties of the cuprates; (ii) on the other hand, Wells et al. H noticed that 
the 2D t — J model prediction fails near X = (0, tt) since the experiments do not confirm 
that the energies of the M and X points are nearly degenerate, as is predicted by t — J model 
calculations ||||. 

It is the purpose of this paper to provide evidence that item (ii) above can be naturally 
accounted for in two different ways, and we shall discuss each of these ideas separately. 
Firstly, if the t — J model is enlarged to contain hole hopping amplitudes at distances larger 
than one lattice spacing, one can indeed fit the ARPES data. There is no symmetry con- 
straining the hole hopping to just nearest-neighbors, and in several studies it has been shown 
that to reproduce the spectra of three-band Hubbard models using one-band t — J models 
on small clusters, it is necessary to include a "t'-hopping" along the plaquette diagonals || . 
As long as the ratio |t'/t| is small, this term does not affect the AFM properties of the model 
near half-filling, since this hopping only connects sites that are on the same sublattice of a 
Neel ordered state. However, such terms will clearly affect the quasiparticle dispersion. 

It has, in fact, become standard to use Hamiltonians including t' hoppings to describe 
several compounds. Some authors J7|] have argued that to fit the ARPES Fermi surface 
of different materials it is necessary to select t' = — 0.2t for La2_ x Sr x Cu04; t' = — 0.45t for 
YBCO; and t' = +0.2t for Nd2_ x Ce x Cu04, and there seems to be some justification for this 
approach. For example, it was recently shown || that the differing magnetic properties of 
the LaSrCuO and NdCeCuO systems could be accounted for only if such a next-nearest- 
neighbour hopping term was included. Since all these compounds have the same Cu0 2 
planes, the strength of these effective "t'-terms" changes among the different cuprates likely 
due to three dimensional effects (bi-layer hopping, influence of apical oxygens, influence of 
orbitals usually neglected in the 3 and 1 band descriptions of the cuprates, etc.). Then, it is 
reasonable to expect that the t — t' — J model can mimic the dispersion of S^CuC^C^ with 
a properly selected amplitude t'. Although it is clear that fitting data with an arbitrary 
parameter is less satisfying than microscopically deducing the value of t', in the absence of 
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a microscopic theory to compute t' the best we can do is show that a reasonable value of 
this amplitude leads to agreement with the experiments. Note that we only have one free 
parameter to fit an entire function of momentum (dispersion), and thus not much extra 
freedom is introduced in the model with a t' amplitude. 
The t — t' — J model Hamiltonian is defined as 

H = t + h.c.) - t' J2 + h.c.) 

<ij> 

+ jECSi-Sj-inmj), (l) 

(ij> 

where c are hole operators, () refers to n.n. sites, (()) refers to n.n.n. sites (along the diagonal 
of the plaquettes), and the rest of the notation is standard. Using the Born approximation 
(extensively discussed in the literature [Hifl) the dispersion of one hole in the t — t' — J 
model was here calculated for several values of t' [[TT|. To select t' the experimentally 
observed isotropy of the dispersion of S^CuC^C^ near the top of the band at (tt/2, ir/2) |§ 
was used. We followed this criterion since it is the most clearly defined feature in the data of 
Ref. || (concentrating our effort on the measured behavior near X would be dangerous since 
for these momenta the ARPES data only provide a bound on the quasiparticle energy). We 
observed that to reproduce the (7r/2,7r/2) isotropy, an amplitude t' ~ — 0.35t is necessary 
(supplemented by J = 0.125eV as energy scale and a ratio J/t around 0.3-0.4, as is usually 
assumed). This value of t' is similar to that found in earlier comparisons based on small 
cluster exact diagonalization results || to the ARPES data [[11]]. The resulting dispersion is 



shown in Fig. 1, and is compared with the experimental data. Also shown is the dispersion 
produced by the t — J model. The agreement near the top of the band is much improved 
by the inclusion of t', especially in the (0, tc) — > (it, 0) direction. It is clear that along this 
direction the t — J and t — t'(= — 0.35t) — J models have substantially different dispersion 
relations. On the other hand, along the (tt, 0) — > (0, 0) line our results have more structure 
than the reported experiments. However, experimentally along this line it is not possible to 
extract a reliable quasiparticle dispersion (see Fig. 2a of Ref. ||) and thus we do not consider 
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this disagreement with states having an energy much less than the chemical potential a 
serious problem. Actually, it is possible that only the quasiparticles near M are relevant for 
normal state and superconducting properties ||12|| . Nevertheless, it is important to fix this 



quantitative discrepancy, and for this purpose we turn to our second technique for producing 
a quasiparticle dispersion that agrees with the data of Ref . . 

The underlying physics of carrier transport in weakly doped CuC>2 planes is governed by 
oxygen hole motion in an AFM background, and thus one might expect that the three-band 
Hubbard model in the strong-coupling limit should be employed. However, as is now 



well known, a renormalization procedure approximately maps the low-energy physics of the 
three-band Hubbard model into that of the one-band t — J model [|14j]. Subsequently, a per- 



turbative calculation [15] prescribed conditions under which further-than-nearest-neighbour 



hopping terms must be included. As demonstrated above, if one wishes to fit the ARPES 
data of Ref. at least the next-nearest-neighbour hopping term t' is required. 

One possible caveat to the reduction to one-band models has been discussed in previous 
papers [|T^-|TB[. To be specific, the reduction from the three-band Hubbard model to the 
t — J model eliminates the spin degree of freedom of the oxygen carriers, effectively producing 
spinless vacancies. However, if the direct oxygen-oxygen hopping is large enough, the oxygen 
hole can acquire a non-zero spin. As a demonstration of the potential inadequacies of any 
renormalization to a one-band model based on carriers possessing no spin, consider firstly 
that unlike the t — t' — J model without the t' term, holes described by the three-band 
model moving in a rigid Ising-like AFM are mobile |17|], and thus these carriers possess a 



spin degree of freedom. However, in order to enhance the derealization of these holes in a 
quantum AFM, the spins become distorted, and the spin of the carrier is quenched |14| , |16 



approximating the spinless vacancies of the one-band t — t' — J model. If one then includes 
a hopping process that is independent of the magnetic background, such as a direct oxygen- 
oxygen hopping term, and the energy scale of this additional hopping term is larger than 
any of the other terms, the spin of the carrier is no longer quenched and one should 
then not expect that such carrier motion (in a strongly correlated magnetic background) 



should be describable by the spinless vacancies of the t — t' — J model. 

Motivated by this observation, we have examined the usefulness of fitting the experimen- 
tal results with the microscopic Hamiltonian of Ref . jH]] . We have studied the analogue of 
the three-band t — J model found from the strong coupling limit of the three-band Hubbard 
model. This is a complicated microscopic Hamiltonian, and we refer the reader to Ref. [T^ 



for a detailed discussion of Eq. (2) of that paper. In addition to this, here we append 
the Hamiltonian of Ref. [ 1(| with a direct oxygen-oxygen hopping term, something that we 



found to be crucial to our success in fitting the ARPES data |§. The parameter values 



for this Hamiltonian jig] can be estimated from the parameters of the three-band Hubbard 



model ||, with the Cu — O exchange interaction reduced by a direct exchange term [19 



Then, following the notation of Ref. |T(| , one finds 

Ji = 0.032, J 2 « 0.16, t a = 0.38, t b = 0.43, 

and t pp = -0.65 (2) 

Extrapolating the analysis of Shastry to include the direct oxygen hopping, and accord- 
ing to the above argument concerning the potential inadequacies of the renormalization to a 
one-band t — t' — J model, the important ratio of hopping parameters is t pp /|t a + tb| ~ 0.8. 
This implies that the renormalization to spinless vacancies may not be entirely justified. 

We have used these parameters in an exact diagonalization study of the Hamiltonian of 
Ref. |Tj| for a 4 x 4 CuC>2 cluster with periodic boundary conditions. Our results are shown 



in Fig. 2. Note that the shape of the quasiparticle dispersion (at least for the small number 
of wave vectors that we can access) is identical to that of Wells' data, which is encouraging. 
However, the bandwidth is roughly half of that found experimentally. This latter result is 
not that surprising given that for a large t pp , important finite-size effects are to be expected 

In order to better assess the quasiparticle dispersion relation that describes a single 
oxygen hole doped into an AFM that follows from Eq. (2) of Ref. [IB"], we have generated 



the quasiparticle dispersion using a semiclassical variational wave function that includes 



the spin fluctuations of the Oxygen holes where the distortions of the Cu spins away from 



AFM order are determined variationally (see Eq. (19) of Ref. (16| for a version of this wave 
function that does not include Oxygen spin fluctuations). We expect that this semiclassical 
representation will lead to an excellent approximation to the true quasiparticle dispersion 
quite simply because the carrier motion has a large direct oxygen hopping. (For example, if 
one modifies the results of Roth |17| to include t pp and allows for the hopping parameters of 



Eq. (2) to become fitting parameters, even if one ignores the Oxygen hole spin fluctuations 
one can again fit the Wells' data (shape as well as bandwidth) reasonably well — this is also 
displayed in Fig. 2.) 



We have included these spin distortions [21,16], and numerically solved the associated 



variational problem to obtain the minimum energy state at each wave vector. This calcula- 
tion has been carried out for an infinite lattice. Here, we have again allowed for the hopping 
parameters of Eq. (2) (above) to be variables that we use to fit the quasiparticle dispersion 
found experimentally, expecting that we obtain values close to those of the perturbation the- 
ory [|IBJ. Our results, for t a ~ 0.0, tb = 0.3, and t pp = —0.35, are shown in Fig. 3. Clearly, 
the agreement with experiment is quite good, and improves that found using the t — t' — J 
model along the (0, 0) — > (0, it) branch. This may point to the importance, at least for these 
crystal momenta, of the spin degree of freedom of the holes. 

Wells et al. || made the apparently reasonable assumption of grouping together the 
results obtained for hole-doped Bi2212 [l|], where flat bands at (0,7r) near the Fermi energy 
were reported, with those for insulating Sr 2 Cu0 2 Cl 2 , concluding that as a function of doping 
states near (0, 7r) must either be created or move up in energy in a non-rigid way. This is one 
possible interpretation. However, the many clear differences in the ARPES Fermi surfaces 
of different high-T c compounds leads us to believe that it is not correct to group Bi2212 
with Sr 2 Cu02Cl2 since these compounds differ in microscopic details that are relevant for 
quantitative comparisons — e.g., the numerical value of the t' amplitude. While it is likely 
that the qualitative mechanism of superconductivity is the same in all hole-doped cuprates, 
their differences in critical temperature and ARPES data must be caused by non-universal 



features. Thus, in our opinion, the data from Bi2212 and S^CuC^C^ cannot be grouped 
together to obtain conclusions about the behavior of correlated electrons as a function of 
doping in the cuprates. To clarify this important point it will be necessary to carry out 
photoemission experiments on the same compound at different densities of carriers. 

Summarizing, using a t — t' — J model it is possible to fit the ARPES quasiparticle 
dispersion data for S^CuC^C^ reasonably well at the top of the valence band. To obtain 
an even better quantitative agreement the use of a strong-coupling three-band Hubbard 
model is required [[U|. 
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Figure Captions 

1. Quasiparticle dispersion of the t — t' — J model calculated in the Born approximation 
for an infinite lattice using t' = — 0.35t (solid line), J/t = 0.3 and J = 0.125eV, com- 
pared against the experimental ARPES data (open circles) for S^CuC^C^ of Ref. [3]. 
The dashed line is the t — J model result of Ref. f5|. 

2. Quasiparticle energies (stars) for an oxygen hole described by the effective Hamiltonian 
of Eq. (2) of Ref. [TO found from an exact diagonalization study on a 4x4 CuC>2 cluster, 



compared against the S^CuC^C^ ARPES data (open circles); the zero of energy of 
the exact diagonalization numbers is offset for clarity. Also, the solid line denotes the 
band of energies for hole motion in a rigid (viz., Ising like) AFM background using 
hopping parameters t a = 0.12, tb = 0.17, and t pp = —0.18. 

3. Quasiparticle energies (solid diamonds) of an oxygen hole in the effective Hamiltonian 
of Eq. (2) of Ref. for an infinite lattice using a semiclassical variational wave 
function that includes the spin fluctuations of the Oxygen carriers. The quasiparticle 
dispersion measured in Ref. [3] is also shown (open cirlces). The hopping parameters 
that produce this quasiparticle dispersion are t a = 0.0, tb = 0.3, and t pp = —0.35. 
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